Quantum Conditions on Dynamics and Control in Open Systems 
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Quantum conditions on the control of dynamics of a system coupled to an environment are 
obtained. Specifically, consider a system initially in a system subspace Ho of dimensionality Mo, 
which evolves to populate system subspaces Hi , H% of dimensionality M\ , M2 . Then there always 
exists an initial state in Ho that does not evolve into H2 if Mo > dMi, where 2 < d < (Mo + Mi + 
M2) 2 is the number of operators in the Kraus representation. Note, significantly, that the maximum 
d can be far smaller than the dimension of the bath. If this condition is not satisfied then dynamics 
from Ho that avoids H2 can only be attained physically under stringent conditions. An example 
from molecular dynamics and spectroscopy, i.e. donor to acceptor energy transfer, is provided. 



The study of open quantum systems [J, 0,0] is of great 
current interest, a consequence of the desire to under- 
stand and manipulate devices that are reliant on quan- 
tum effects. Of particular interest is the control of such 
quantum systems, with coherent control 0, [1| being a 
most promising approach. 

The essential principle of the coherent control of 
atomic and molecular processes relies on the creation 
of multiple interfering pathways to the same final state 
0, @| • Manipulation of laboratory parameters that af- 
fect these coherent pathways then allows direct control 
over the associated interference contributions, and hence 
control over the relative cross sections involved in the 
process. Numerous theoretical scenarios and experimen- 
tal results, in a host of applications, have successfully 
demonstrated such control. 

Despite the importance and generality of the field, 
there are remarkably few formal results that provide nec- 
essary and sufficient conditions under which complete 
control of the system dynamics is possible. These in- 
clude theorems such as that of Huang- Tarn-Clark @, a 
theorem by Ramakrisnan et al. on the dimensionality 
of the Lie Algebra induced by the interaction between 
the system and the control field 0], and a theorem by 
Shapiro and Brumer Q, where control was shown to de- 
pend on the dimensionality of the controlled subspaces. 
These theorems all deal with closed systems, i.e. those 
that are not in contact with an environment. Real quan- 
tum systems are, however, not isolated. They often inter- 
act with a surrounding environment, which is regarded 
as an uncontrolled system of arbitrarily large dimension. 
As a result the system is open to its environment, and 
proofs regarding control in open systems are becoming 
the focus of serious attention. For example, Reference 
Q provides formal theorems on controllability in open 
system dynamics treated via the Kraus representation. 

In this letter we significantly generalize a prior 
argument Q for closed systems and obtain quantum con- 
trol conditions in open systems. We show that the con- 
dition for a specific control problem in open systems is 
stronger than that in closed systems, but that the rela- 



tive dimensionality of the different subspaces remains a 
crucial determinant of controllability. In addition, our 
approach gives the solution to the control problem when 
it is achievable, and provides a constructive method for 
extending various closed system proofs to the open sys- 
tem case. 

To expose the difference between the open and closed 
cases, to adopt a common notation, and to provide 
a generalized formalism, we first reconsider the closed 
system 0. Consider a system whose Hamiltonian eigen- 
states are partitioned into bases for three subspaces: Ho, 
Hi, H 2 , of dimensionality M , Mi and M 2 , respectively. 
Let the basis vectors \i, m) span the subspace Hi, where 
i = 0,1,2 labels the subspace and n; = 1, ...,M, labels 
the states in Hi. The total dimension of the system is 
M = Mq + Mi + M-2- A superposition of basis vectors 
from Ho is chosen as the initial state. We assume that 
the system initially resides in Mq and that the state can 
flow into the other two subspaces under the dynamics. 
The question under consideration is: under what condi- 
tions can one prevent dynamics into the subspace H 2 , 
by preparing states in Hq that go solely into H{1 To 
resolve this problem, Ref. [H considered linear superpo- 
sitions of states in Hq and utilized interference between 
the resultant paths to Hi and H 2 in order to control the 
dynamics. 

Using the evolution operator, U we have: 



W)=U\i> ) 



(1) 



where \ipo) is the initial superposition in Mq and \^p) de- 
notes the final state. With U represented as an M x M 
dimensional matrix, and \ip) expressed as an M column 
matrix we can examine the number of simultaneous equa- 
tions that need to be solved for \tp) to have zero popu- 
lation in the rows corresponding to H 2 . Using this ap- 
proach, Ref. Q obtained conditions on the relative di- 
mensionality of the different spaces by requiring 



(2,k 2 \U\i> ) = 



(2) 



for all states k 2 in H 2 . Specifically, they showed that 
if Mq > M 2 , it is always possible to prevent transi- 
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tions from Hq into H 2 ■ By contrast, a large number 
of difficult-to-satisfy linear dependence conditions are re- 
quired to prevent transitions if M ^ M 2 . Significantly, 
the result only depends on the dimensionality M and 
M 2 of these two subspaces. 

Consider then the same system, but where it is now 
open to an environment, which can be either finite or in- 
finite dimensional, and that may, or may not, be in the 
presence of external fields. As is typically the case, we as- 
sume that the system can be addressed (e.g. via laser ex- 
citation) to prepare the initial state, but that the environ- 
ment remains unaffected. Further, we assume that the 
initial system-bath density matrix ptot(O) = p(0) <8>ps(0) 
is a separable product of the initial system p(0), usually 
a pure state that can be selectively prepared, and initial 
bath ps(0) density matrix[l(|. The open system after 
evolution can be characterized by an M x M reduced 
density matrix p, defined as p =TvBPtot, where ptot is 
the final time total system-bath density matrix and Tr^ 
indicates a trace over the bath. Focusing on the reduced 
density matrix (i.e. ptot traced over the bath) allows us to 
continue to use the concept of dimensionality in an open 
system. In general, the reduced density matrix evolves 
according to a linear transformation 



p = Y J A a P {Q)B P 

a, /3 



(3) 



For instance, in the case of natural quantum evolution, 
denoting the propagator for the full (system + bath) as 
U, we have that 



Ptot = Up tot (0)tf = Up(0) ® PB (0)U\ 



(4) 



where we have used the separability of ptot at time zero. 
In general Pb(0) — J2 a bP ab \ a ) (^l> where |o) and \b), and 
the e) below, are eigenstates of the bath Hamiltonian. 
By tracing Eq. (j4|) over the bath coordinates, we have 
that 



P : 



53 A a ' e p{0)B 



a.b.e 



(5) 



with A a - e = (e\U\a) and B a ^ e = p ab {b\W\e). 

Under the assumption of complete positivity ll| , Eq. 
(|3|) simplifies and the reduced density matrix evolves 
according to the operator-sum representation (i.e. the 
canonical Kraus representation), which is generic physi- 
cally for the initial product density matrix [3j]: 



= J2 £"/°(0)£ Qt 



(6) 



Here the sum over Kraus operators [Tl| E a is such that 
Ei=i E af E a < I. Specifically, £l=i E^E<* = I 
for trace-preserving quantum operations and 
Ea=i E a E a ^ < I for non-trace-preserving quan- 
tum operations such as quantum measurements. Note, 



significantly, that the operator-sum representation has 
at most M 2 Kraus operators, i.e., 2 < d < M 2 and that 
the operators are fixed by the given system, bath and 
any incident external fields. 

We continue below to consider the general form [Eq. 
©], but subsequently focus on evolution under the Kraus 
representation. 

Note first that Eq. ([3]) means that the operator p(0) 
can be linearly transformed into the operator p. If we 
denote the operation connecting the matrix p(0) to p by 
V, then V, in a particular representation, is a matrix with 
four subscripts. That is, 



Pts 



53 Vtst's'Pt' 



'(0) 



where t or s denotes the indices [i, rii) of the basis vectors 
\i, rii) . Rewriting p as M 2 dimensional column vector de- 
noted p, Eq. ([6|) can be equivalently rewritten as: 

P=Vp(0) (7) 

where the explicit form of the M 2 x M 2 matrix V is 

v = 53 A ° ® B& 

Equation jJJ) is seen to have the same form as Eq. |lj , 
allowing the approach used earlier for a closed system to 
be extended to open systems. Consider then an initial 
reduced density matrix that is in Hq, i.e. 

/ Pul 
p(0) = 
V 

where p M 2 is an Mq x Mq matrix. In order to prevent 
dynamics from going into H 2 , we require 



P(2,fe 2 )(2,fe 2 ) 







(8) 



For the case of interest the Kraus representation ap- 
plies so that a propagated initial pure state p M 2 = 
IV'o) (V*o| j becomes 

d 

P(2 M )( m) = 53 (2, k 2 \ E a |^o> (M E^ a \j, n 3 ) 

Q = l 

Given that the diagonal element P(2,k 2 )(2,k 2 ) ^ s 

E«=i I <2, fc 2 | £7" I -0o > 1 2 , satisfying p (2 ,fc 2 )(2,fe 2 ) = 
means requiring 



(2,k 2 \E a \i; ) =0 



(9) 



This being the case, demanding zero population in H 2 
also implies that all elements P( 2 ,k 2 )(j.n j )j which include 
coherences with the H 2 Hilbert subspace, are also zero. 

Equation © has a nontrivial solution if Mo > dM 2 , 
and control via initial state preparation is achievable no 
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matter what the form or dynamics of the ci-dimensional 
Kraus operators. Further, this equation provides the ini- 
tial state that allows for the desired control. However, 
this control condition is, as expected, far more stringent 
in the open system than in the closed system case, where 
d = 1. Note, however, that the bath effects are still lim- 
ited, The minimum d can be[l2| as small as two, and the 
maximum d is M 2 , which can be far smaller than the 
dimensionality of the bath. 

If Mq < dM2, control via initial state preparation is 
also far more difficult in the open system than in the (al- 
ready difficult) closed system case. Define a dM 2 x Mq di- 
mensional matrix W with matrix elements W '12k, a), (On) — 
(2, k\ E a |0, n) . In the case where Mq s$ dM 2 , the rank of 
W is equal to Mq unless all Mq x Mq dimensional subma- 
trices of W are singular. Hence, non-trivial solutions to 
Eq. (|9|) exist if det(w|^ M ) = 0, where (A;) numbers all 
of the submatrices. This condition also implies that a set 
of columns of W are linearly dependent. Hence control is 
possible for a specific class of W, expected to be difficult 
to obtain physically, 

Note that the definition of the subspaces Hi and as- 
sociated dimensions Mi can, in some instances, also be 
manipulated. For example, in bound state systems (such 
as that in the example below) these subspaces can be de- 
fined by accessing only specific system eigenstates using 
a pulsed laser field. This facility might prove additionally 
useful in attempting to satisfy the Mq > dM 2 bound of 
this theorem. 

These results allow numerous applications. Consider, 
for example, electronic energy transfer from donor to 
acceptor molecules, where both may be part of one 
larger molecule. These systems are ubiquitous, rang- 
ing from relatively small systems 13] to large structures 
such as carotenoid-to-bacteriochlorophyll energy transfer 
in photosynthesis 14( . The most interesting cases take 
place in condensed matter environments, i.e. open sys- 
tems. Studies of the dynamics and spectroscopy of such 
systems can be carried out by laser excitation of the sys- 
tem from a lower electronic manifold (here H ) to the 
donor (here Hi). Electronic energy transfer from the 
donor (Hi) to the acceptor (H 2 ) is then measured. Here, 
the subspaces Hi, and associated dimensions Mi, are de- 
termined by a combination of the molecular state den- 
sities and the width of the laser pulses that prepare a 
preliminary superposition of states on Hq and that sub- 
sequently excite the system into Hi and H 2 . 

In some systems of interest, excitation from Hq to the 
donor is contaminated by partial excitation of the accep- 
tor as well, with a concomitant reduction in the quality of 
the data on the subsequent electronic energy transfer dy- 
namics. A considerable improvement would result from 
being able to excite Hi with reduced population trans- 
fer to H 2 from H . Results of the open system theorem 
indicate that this is (a) difficult to achieve physically if 
M < dM 2 , and (b) attainable if M > dM 2 . In the 



latter case, one could ensure significantly reduced, and 
ideally zero, acceptor population at the target time. The 
extent to which this is achievable is dependent upon the 
particular system; specific systems of this type will be 
the subject of future study. 

In summary, the quantum conditions obtained above 
are completely general, applying to both systems that 
are controlled, as well as to uncontrolled system evo- 
lution. The result establishes an important inequality 
between the dimension of the subspace M 2 which we de- 
sire not to populate, the initial subspace Mq from which 
dynamics evolves, and the dimensionality of the Kraus 
representation. It is a dynamics-independent property of 
rf-dimcnsional Kraus evolution. As long as the geometric 
condition Mq > dM 2 is satisfied, control via initial state 
preparation is achievable. Further, the fact that control 
can be achieved via initial state preparation, in the pres- 
ence of an environment under the prescribed conditions, 
is useful for control applications in realistic systems. 

Two supplementary remarks are in order. First, 
this theorem places emphasis on the importance of the 
dimensionality d of the Kraus representation, known 
to satisfy 2 < d < [Mq + M 1 + M 2 ) 2 . Hence, these 
results should motivate further studies to determine the 
actual d for realistic systems. Second, we note that an 
analogous method to that described above can be used 
to extend other closed system proofs, such as those of 
Refs. @ and 0, to open systems. 
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